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Abstract 

We consider a perturbation of the Anosov-type system, which leads to the 
appearance of a hierarchical set of islands-around-islands. We demonstrate by 
simulation that the boundaries of the islands are sticky to trajectories. This 
phenomenon leads to the distribution of Poincare recurrences with power-like 
tails in contrast to the exponential distribution in the Anosov-type systems. 
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1. INTRODUCTION 

Anosov maps (or Anosov-type maps, or Anosov diffeomorphisms) can be consid- 
ered as a strong idealization of realistic systems. The reason of this is that Hamilto- 
nian systems are not uniformly hyperbolic and their phase space features a complex 
mixture of zones of different dynamics. The boundaries of the zones influence the 
trajectories and drastically change asymptotic properties of transport. The ways how 
the Anosov systems can lose their hyperbolicity is the subject of various investigations 
[H El El El E] The main questions discussed here are a description of the topological 
pattern of the destroyed hyperbolicity and the description of the consequences of that 
loss. For the latter issue, it could be, for example, the changes in particles transport 
due to the changes in the phase space topology that occurs in Hamiltonian systems 
with mixed phase space jHj, i.e. phase space with chaotic sea and islands. 

In this paper we make few comments related to the destruction of hyperbolicity 
and to the problem of transport. The considered system is similar to the introduced 
one in pp. It is a map that acts on the torus T 2 and it depends on an intrinsic param- 
eter K G (—00, 00), and on the perturbation parameter €q. The uniform hyperbolicity 
exists when eo = and in a region K + of the values of K. The system is most sensi- 
tive to the perturbation in the vicinity of edges of K + . In the extended space T 2 x K 
the system is unstable with respect to an arbitrary small eo- Increasing eo leads to a 
sequence of bifurcations. We demonstrate by simulations the appearance of a hierar- 
chical set of islands in the stochastic sea and a stickiness of their borders. In other 
words, trajectories spend enormously long time in the vicinity of the islands' border 
making the transport to be anomalous. The "anomalous" means that the wandering 
of trajectories strongly differs from the random walk in the case of uniform hyperbol- 
icity. This property is established numerically by applying the method of e-separation 
along reference trajectories developed in [HEH^]. The mam result can be expressed 
in a simplified form as follows: After the loss of hyperbolicity, the considered Anosov 
system reveals the power law of the distribution of Poincare recurrences when t (time) 
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— > oo, in contrast to the exponential distribution for unperturbed Anosov system |lUj . 
2. DESCRIPTION OF THE SYSTEM 

Consider a dynamical system (x, y) e T 2 defined on the two-dimensional torus by 
the equation of motion 

Vn+i = y n + f(x n ] eo, K); x n+1 = x n + y n+1 , (modd 1) (1) 

and a differentiate function 



f(x n ; e , K) = Kx n - e ^(x n ) 



(2) 



with two parameters K and eo- The map (0) is area-preserving for any f(x;eo,K). 
While the map acts in the phase space which is T 2 , it doesn't transform the boundary 
of the torus dT 2 into dT 2 . Particularly, for eo = it is evident since K is not integer. 
This is a typical physical situation. One can interpret the map (0) as a generator 
of trajectories {xo, yo; X\, y\\ ...} in the lifted space (x, y) G (— oo, oo; — oo, oo), which 
should be, afterward, bended and put inside the compact phase space (x, y) moddl. 
This way of the consideration of the map (0) permits to study infinite diffusion in 
the lifted space and local measures, Poincare recurrences measure, etc. in the phase 
space T 2 = (x,y) moddl. 

The eigenvalues of the tangent to (JTJ) map are 



Al,2 : 

For eo = the region 



l + \K- \e^'{x) ) ± 



l+ l -K- l -e^'{x)^ -I 



1/2 



(3) 



K + : K < -4 , K > ; (e = 0) 



(4) 



corresponds to uniformly hyperbolic dynamics, and the region 



K : -4 < K < 



(5) 



corresponds to periodic dynamics with zero Lyapunov exponent. For our following 
consideration a specific form of ip(x) is not so important. Nevertheless, we specify it 



for a convenience 

ijj(x) = smx ■ cos 2 x , (6) 

as it was considered in [H Hj . 

The system (JTJ), © represents Anosov diffeomorphism for e = 0, K G K + , and 
it is structurally (topologically) stable for e — > 0. Our goal is to consider structural 
stability in the extended space T 2 x K where K G K. 
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Figure 1: Poincare map of a trajectory for the system (1) with (2), (6) and K = 0.2, 
e = 0.2 (left), e = 0.16 (right). 

3. ISLANDS AND ISLANDS HIERARCHY 



Property 1. For arbitrary small eo there exists Kf G K such that for K G 
(K~, K+) phase space T 2 has islands of periodic orbits that appear as a result of the 
bifurcation. The main reason for this is that the hyperbolicity has edges defined in 
(jU, ^and arbitrary small eo influence the system's stability fairly close to the edges. 
An example of the island born is in Fig. ^ For eo 7^ the domain of stability 



K 



-4 + e ip'(x) < K < e V>'(z) 



(7) 



depends on x, and the phase space becomes non-uniform. 

For small e the most interesting domain is near x = and small K. For small x 
we have from (J2J): 

fix) = 5x + (7/6)e x 3 + 0(x 3 ) ; 5 = K - e . (8) 

The map ( (JTJ), in the vicinity of (0,0), gives for 5 — > 

7 

y n+2 = (1 + 5)y n + 25x n + -e x 3 n + 0(5 2 ) 

o 




that is area-preserving up to the terms of the order <5 3//2 , Sx. On the same level of 
accuracy, the map (JHJ) can be written as Hamiltonian equations 



V = 



dH 

dx 



x = 



dH 

dy 




n 



(10) 



These expressions can be easily improved to have terms of higher order but we will 
not need it. An elliptic island appears around the point x = y = for 5 < 0, i.e. for 
K < e in correspondence to (J7J) and ©. The island has a scaling 



The value 5 = is the point of degeneracy and the bifurcation "opens" an island 
in the phase space [I]. Increasing of \8\ leads to a sequence of bifurcations that change 
the phase space topology. A particular case, shown in Fig. EJ is of special interest. 
For the specific values of the pair (K, e) we have a hierarchical structure of islands- 
around-islands with a sequence of islands: 1—4— 4-3— 4-3-8— 4-3-8-8— 4-3-8-8-8 — . . .. 
The values 



are found numerically, for which the hierarchy of islands has a constant asymptotic 
proliferation number 8. Sequences of islands with bounded values of rrij : rrij G 
(m~,m + ) is called a hierarchy of islands since the full number of islands at the 
hierarchy level n is equal m" with some rh e (m~, m + ) and n — > oo. 

Property 2. There exists an infinite number of hierarchy of islands sequences 
(1, mi, m 2 , . . .) and the corresponding pairs (K, eo). For a fixed e (or K), it was 
assumed in [TT] (see also jHj), that the values of K (or eo) to have an island hierarchy, 



y~\5\ l ' 2 x, 



Of-o) 



(ii) 



(K*,e*) : K* — I 



1.689538... 



(12) 
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Figure 2: Sticky islands for the same equation as in Fig. ^ with K = K* = 1 and 
eo = £q = 1.689538. (a) full phase space; (b) zoom of the group of islands on the top 
right of (a); (c) zoom of a next generation island from the left island of (b); (d) zoom 
of the lower island in (c). 
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are at least as dense as rationals. Many examples of this type of islands topology can 
be found in the review 

4. STICKINESS AND QUASI-TRAPS 

The phenomenon of a long time stay of particles in the vicinity of the islands border 
is called stickiness. While it is known that the stickiness is due to the presence of 
cantori and some renormalization group approach can be applied to study particles 
transport through the cantori IT?] , a rigorous theory of stickiness does not exist yet 
and many features of the stickiness are waiting for our understanding (see more in [S|). 
A definition of stickiness can be provided using the notion of dynamical traps^3], or 
for Hamiltonian systems, quasi-traps. 

Let T be the compact phase space of a dynamical system that is Hamiltonian and 
chaotic in a part of T, and let C T be an island imbedded in the stochastic sea. 
Consider a border of the island, OF a, and an annulus A = OTa x L with a fairly 
small width i. For any initial condition (x, y) G A, one can consider a time t(x, y; A) 
of a trajectory first exit from A. This time is non-zero and non-infinite, since the 
dynamics is area-preserving. 

The function p(r; A) is called a distribution function of the trapping or escape 
time from the domain A if 

/•oo rt 

V{t;A)= drp(r; A) = 1 - / drp(r; A) (13) 
Jt Jo 

is the probability to escape from A after time > t. Consider moments 

(r q ) A = / drr*p(f;A) . (14) 
Jo 

For q = 1 we have (t)a < oo due to the Kac lemma |14j . i.e. the mean trapping time 
is finite. The domain A is called quasi-traps if there exists such q > 1 that {t q )a = oo 
[iniini- Evidently the distribution function pit; A) should satisfy the condition 

p(t;A)~l/ t \ 7 >1, (t-oo), (15) 

for quasi-traps. 
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There are many examples of systems with quasi-traps: standard map, Sinai and 
Bunimovich billiards (with 7 w 3), rhombic billiard (with, most probably, 7 close to 
2), and others (see in review [Hj). It was assumed in jHj that realistic Hamiltonian 
systems with a hierarchy of islands reveals the dynamics with quasi-traps. In this 
paper we demonstrate by simulations that similar property exists for the perturbed 
Anosov diffeomorphism (JTj), ((2j) near the point (|T2|) where a hierarchy of islands 
appears. 

Visualization of a trajectory of the map (1) in Fig.|2]shows stickiness as dark strips 
around the islands' borders of all presented generations. Although the stickiness 
is similar to one, observed in [llj . it is not easy to follow the same methodology 
since strong mixing for the Anosov-type system. In order to provide a quantitative 
characteristic of the stickiness, we apply a method of e-separation of trajectories 
developed in jZUHlEl- We are starting a pair of trajectories in some domain B C T 
fairly far from the islands' area A. One trajectory of the pair will be called the 
reference trajectory, and the initial distance between trajectories is d <C dp, d r is a 
diameter of T. The distance at time t is d(t). Let t\ be a first instant when separation 
of the trajectories reaches the value e that satisfies the condition 

d < e < d r . (16) 

After that, at time instant ti we neglect the trajectory that is not the reference one 
and start another trajectory at a distance do from the reference one. This new pair 
will be e-separated first time at t 2 > t\. Such a procedure continues as long as the 
reference trajectory runs. 

As a result we collect different intervals {r,} = (ti,t 2 ,...), tj = tj+\ — tj of 
e-separation along the reference trajectory. The difference between this procedure 
and the simulation of Lyapunov exponents is two-fold: it is in the condition (16), 
since e -C dr, and the ratio e/d Q is fixed. Finally, from a large number of the initial 
reference trajectories in B, we obtain a set {t}n = {ti, ■ ■ ■ , Tat} that provides a 
frequency function Pat(t; B, A) of the relative number of events when r e (r, r + Ar). 
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If the number of elements in {r} is fairly large, then Ptv(t; B, A) does not depend on 
B, providing Tb H Ta = 0, and it reaches some limit 

p(r; A) = lim p(r; A, B) = lim N(r; A, B)/N , / dr p(r; A) = 1 (17) 

N^oo iV^oo JO 

where N(t; A, B) is the full number of e-separations with separation time r G (r, r + 
rfr) from the set {t}n, and initial conditions (xo,yo) £ B. 

A typical distribution N(n), where n is discrete time and A, B are omitted, 
is shown in Fig. |3] at the degeneracy point K = e - F° r the initial distance in 
pairs we select do = 10~ 6 and e = 10~ 2 . All 4048 reference trajectories were run 
for 10 7 iterations and full number of the separation events iV ~ 10 9 . There are 3 
different regions in the behavior of N(n). Before the sharp peak at n p ~ 100 there 
is N(n) ~ exp(ara), where a is of the order of the mean Lyapunov exponent (see 
more below). After the sharp peak, within the interval log 10 n G (2,2.2) there is an 
exponential decay of N(n), and for log 10 n > 2.3 the function N(n) is close to the 
polynomial decay. The latter statement is provisional due to the lack of statistics 
and small time interval. Because of inhomogeneity of phase space, the finite time 
Lyapunov exponent at has a distribution function, and the first two domains with 
the peak fit the corresponding theory of the probabilistic analysis of the distribution 
of a t (see for example in^l)- The third domain shows a new feature that can be 
interpreted as influence of the bifurcation point that leads to the appearance of the 
stickiness near (x,y) = (0,0). 

Another situation occurs for the values of parameters (K*, e*), given in (|T2~jl when 
the stickiness to the hierarchical islands is strong. For this specific case, let us denote 

Pin) = lim N(t, A,B)/N (18) 

(compare to ()17j)). A corresponding result is shown in Fig. 0] Within four decades of 
n there is a clear power- wise behavior of P(n): 

P[n) ~ const. n~" f , (n — > oo) (19) 
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with 7 f=a 3.3. To analyze this result, we need to 

Conjecture j^]: Assume that the phase space F of compact dynamics can be split 
in two parts: A that corresponds to the only sticky (singular) zone and T\A, and 
that A R T\A = 0. Then for fairly small B C T\A the following equality is valid 

Tree = 7 > (t ^ OO) (20) 

where 7 rec is the recurrence exponent of the distribution of Poincare recurrences to B 

-Prec(^) — const. n _7rcc 

oo 

Prec(n) = 1 . {n oo) (21) 

n=l 

In other words, distribution of e-dispersion of trajectories for small e is defined asymp- 
totically by the sticky zone. 

The value 7 rec = 3.3 may be of a universal character since similar value was 
obtained for the standard map in [11]. It is worthwhile to check the statements (|20j) 
and (|21|) straightforward by calculating the distribution of Poincare recurrences. The 
corresponding data are given in Fig. For the time interval n < 10 4 

•Prec(rc) ~ exp(-n/n) (22) 
but for n > 10 4 there is a crossover to the power law (}2T]) with the same exponent 

7rec = 3.3. 
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Figure 5: Distribution of Poincare recurrence time for special values K* — 1, — 
1.689538 
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5. CONCLUSIONS 

Our numerical results demonstrate that the perturbed Anosov map shows the 
features of the typical realistic Hamiltonian systems of low dimension: appearance 
of islands as a result of a bifurcation; sequence of topological changes in the phase 
space when we change parameters of the system; islands-around-islands hierarchical 
set; stickiness of the islands' boundaries. The most important new features of the 
perturbed Anosov map is the transform of the Poincare recurrences distribution from 
the exponential decay to the power- wise decay as t —>■ oo. Such a change is important 
for numerous applications being responsible for different anomalous properties of the 
observable macroscopic characteristics [6j. 

Both, Property 1 and Property 2 formulated in Sec. 3, include as their main 
part the existence of different infinite hierarchies of islands. These statements are 
conjectures based on high accuracy simulations. The area of stickiness of orbits can 
be characterized by a low rate of dispersion of trajectories. That means that the long 
time, necessary for trajectories separation, can be effectively lower if the accuracy 
is insufficient. A possibility to observe a crossover from exponential dispersion of 
trajectories to the power one just confirms the high level of accuracy of simulations. 
To be more specific, let us mention that the critical value of perturbation is defined 
up to six decimal digits in (|12j). The crossover from the exponential decay of P(n) to 
the power one occurs at log 10 n « 4 (see Fig. EJ) and the power behavior of P(n) is 
up to log 10 n > 6. These results are consistent with level of accuracy of e^. 
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